Quintessence models of Dark Energy with non-minimal coupling 
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We explore quintessence models of dark energy which exhibit non-minimal coupling between the 
dark matter and the dark energy components of the cosmic fluid. The kind of coupling chosen is 
inspired in scalar-tensor theories of gravity. We impose a suitable dynamics of the expansion allowing 
to derive exact Friedmann-Robertson- Walker solutions once the coupling function is given as input. 
Coupling functions that lead to self-interaction potentials of the single and double exponential types 
are the target of the present investigation. The stability and existence of the solutions is discussed 
in some detail. The models are tested against the observational evidence. 

PACS numbers: 04.20. Jb, 04.20.Dw, 98.8Q.-k, 98.80.Es, 95.30.Sf, 95.35.+d 



I. INTRODUCTION 

At the time, all the available observational evidence 
suggest that the energy density of the universe today is 
dominated by a component with negative pressure of, al- 
most, the same absolute value as its energy density. This 
mysterious component, that violates the strong energy 
condition (SEC), drives a present accelerated stage of 
the cosmic evolution and is called "dark energy" (DE). 
A pleiad of models have been investigated to account for 
this SEC violating source of gravity. Among them, one 
of the most successful models is a slowly varying scalar 
field, called 'quintessence' [1,2]. 

Many models of quintessence assume that the back- 
ground and the dark energy evolve independently, so 
their natural generalization are models with non-minimal 
coupling between both components. Although experi- 
mental tests in the solar system impose severe restric- 
tions on the possibility of non-minimal coupling between 
the DE and ordinary matter fluids [3], due to the un- 
known nature of the dark matter (DM) as part of the 
background, it is possible to have additional (non gravi- 
tational) interactions between the DE and the DM com- 
ponents, without conflict with the experimental data.* In 
the present paper we will be concerned, precisely, with 
non-minimal coupling between the DM and DE compo- 
nents of the cosmological fluid. Since, models with non- 
minimal coupling imply interaction (exchange of energy) 
among the DM and the DE, these models provide new 
qualitative features for the coincidence problem [5,6]. It 
has been shown, in particular, that a suitable coupling, 
can produce scaling solutions. The way in which the 
coupling is approached is not unique. In reference [5], 
for instance, the coupling is introduced by hand. In [6] 



Tt should be pointed out, however, that when the stability 
of dark energy potentials in quintessence models is considered, 
the coupling dark matter-dark energy is troublesome [4]. 



the type of coupling is not specified from the beginning. 
Instead, the form of the interaction term is fixed by the 
requirement that the ratio of the energy densities of DM 
and quintessence be a constant. In [7], a suitable inter- 
action between the quintessence field and DM leads to 
a transition from the domination matter era to an ac- 
celerated expansion epoch in the model of Ref . [6] . The 
authors obtained new features for the coincidence prob- 
lem. A model derived from the Dilaton is studied in [8]. 
In this model the coupling function is chosen as a Fourier 
expansion around some minimum of the scalar field. 

A variety of self-interaction potentials have been stud- 
ied. Among them, a single exponential is the simplest 
case. This type of potential leads to two possible late- 
time attractors in the presence of a barotropic fluid [9,10] : 

i) a scaling regime where the scalar field mimics the dy- 
namics of the background fluid, i. e., the ratio between 
both DM and quintessence energy densities is a constant 
and ii) an attractor solution dominated by the scalar 
field. Given that single exponential potentials can lead to 
one of the above scaling solutions, then it should follow 
that a combination of exponentials should allow for a sce- 
nario where the universe can evolve through a radiation- 
matter regime (attractor i)) and, at some recent epoch, 
evolve into the scalar field dominated regime (attractor 

ii) ). Models with single and double exponential poten- 
tials has been studied also in references [11,12]. 

The aim of the present paper is to investigate mod- 
els with non-minimal coupling among the components of 
the cosmic fluid, in which, single and double exponen- 
tial self interaction potentials are involved. To specify 
the coupling, we take as a Lagrangian model, a scalar- 
tensor theory with Lagrangian written in Einstein frame 
variables. In this paper, additionally, we impose the dy- 
namics of the expansion by exploring a linear relationship 
between the Hubble parameter and the time derivative of 
the scalar field. Using this relationship we can solve the 
field equations explicitly. A flat Friedmann-Robertson- 
Walker (FRW) universe filled with a mixture of two in- 
teracting fluids: DM and quintessence, is considered. 

The paper has been organized as follow. In section 2 
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the details of the model are given. The method used to 
derive FRW solutions is explained in section 3. We use 
the approach of 'two fluids': a background fluid of DM 
and a self-interacting scalar field (quintessence). Two 
different couplings that lead to self-interaction potentials 
of single exponential and double exponential class are 
studied separately. In section 4 a study of existence and 
stability of the solutions is presented. The models studied 
in the former sections are tested against the observational 
evidence in section 5. Finally, in section 6, conclusions 
are given. 



II. THE MODEL 



III. DERIVING SOLUTIONS 



Summing up equations (3) and (4) one obtains: 
H + 3H 2 = ^—^ Pm + V. 



(8) 



Let us fix the dynamics of the cosmic evolution by intro- 
ducing the following constrain 



which implies 



a = e*/\ 



(9) 



(10) 



We consider the following action inspired in a scalar- 
tensor theory written in Einstein frame, where the 
quintessence (scalar) field is coupled to the matter de- 
grees of freedom: 



S = J m d^^\{ l -R- l -{V^f 

-V^) + C 2 (^)£ (matter) }. 



(1) 



In this equation R is the curvature scalar, C 2 {<p) - the 
coupling function, and C ma tter is the Lagrangian density 
for the matter degrees of freedom. 

We are concerned here with a flat FRW universe that 
is filled with a background prcssurclcss dark matter and 
a quintessence field (the scalar field (f>) . The line element 
is: 



ds 2 



-dt 2 + a(t) 2 6 il dx i dx j , 



(2) 



where i,j = 1,2,3 and we use the system of units in 
which 8irG = c = Ti = 1. The field equations that are 
derived from the action (1) are: 



3H 2 = p, 
2H + 3H 2 = (l-7)p. 



< + 3H4>= -V - {\nX)'p n 



(3) 

(4) 
(5) 



where we have introduced the reduced notation X{(f) = 
C(0)(37-4)/2 Thc " continuity " equations for back- 
ground (DM) and the quintessence field are: 



p m + 3~fHp m = -(]nXy<j> p m , 
or, after integration 

Pm = Ma-^X-\ 



(6) 



(7) 



where M is a constant of integration. In the former equa- 
tions thc dot accounts for derivative in respect to the 
comoving time t, while the comma denotes derivative in 
respect to (f>. 



and rewriting of (7) in the form: 

p m {(j>) = Me-T*!- 1 ^). 



(11) 



If one substitutes Eqn.(9) in (5) and the resulting equa- 
tion is compared with (8) then, one obtains a differential 
equation on the potential V, that relates V and the cou- 
pling function X (and their derivatives): 

dV „ r . „ , fdlnX A( 7 -2)\ 

If one realizes that d<j) = \d(\na) (see equation (10)), 
then the equation (12) can be rewritten as: 



x 



(13) 



where now the comma denotes derivative in respect to the 
variable N = In a. With thc constrain (9) and the equa- 
tion (11) we can integrate the equation (3) in quadra- 
tures: 



d<t> 



2A 2 



yjMe-^X-^+Vi^) V 6 A ' 



(t + *o), (14) 



or, if one introduces the time variable dr — 
e -^k<$>X- x l 2 dt, (14) can be written as: 



/ 



2X 2 



6 - A 2 



(r + r ). (15) 



In consequence, once the function X((f>) (or X(a)) is 
given as input, then one can solve equation (12) (or (13)) 
to find the functional form of the potential V(</>) (or 
V(a)). Thc integral (14) (or (15)) can then be taken 
explicitly to obtain t = t(<fr) (or r = t (</>)). By inversion 
we can obtain <f> = <f>(t) (or <fi = (j>{r)) so, the scale factor 
can be given as function of either the cosmic time t or 
the time variable r through Eqn.(10). 
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A. Particular Cases 

In this section we study separately two different choices 
of the coupling function X(<f>) that lead to double and sin- 
gle exponential potentials respectively. The importance 
of this class of potentials in cosmology has been already 
outlined in the introductory part of this paper. 



1. Case A 



Let us consider the simplest case when 

4r (InXO)) = const = a=> XU) = A e Q ^, (16) 
d(f> 

where Xq is an integration constant. In this case, by 
integration of (12) one obtains 



V{4>) = Voe-^ + Woe-^ 3 -*/^, 

where the constant 

2M (2a- A(2- 7 ) s 



Wo 



X \a + 37/A - A 



(17) 



(18) 



We are faced with a self-interaction potential that is a 
combination of exponentials with different constants in 
the exponent. The usefulness of this kind of potential 
has been already explained in the introduction and will 
be briefly discussed later within the frame of the present 
model, when we study the stability of the corresponding 
solution. If one introduce the variable z = e~2^ 7 the 
equation (15) can be written as 



/ 



dz 



Vz 2 + a 2 



2X 2 X Q V 

— ( T + To ) (19) 



A 2 



a" 



where a 2 = (X V )/(M + W X ) and I = a + 3 7 /A - A. 
Integration of (19) yields to: 



<Mt) = + Zra(sinh[ / u(r + T )]) 2/1 
and, consequently: 

a(r) = a (sinh \jj,(t + T )]y 2/l . 



(20) 



(21) 



where fi = -j\J -0zi^r-- On the other hand, by us- 
ing the Fricdmann equation (3) and inserting therein the 
potential (17), one obtains the Hubble parameter as func- 
tion of the quintessence field <fi: 



H 2 W 



2V a 
6 - A 2 



b 2 e 



-(a+3 7 /A)0 



-\<t> 



(22) 



where we have considered Eqn.(ll). In terms of the scale 
factor one has, instead, 



H 2 (a) 



2V 
6- A 2 



fe 2 a -(aA+3 7 ) +a -A 2 



(23) 



The former magnitudes can be given in terms of the red- 
shift variable if one considers that a(z) = 00/(1 + 2), 
where ao = a(z = 0) (for simplicity of the calculations 
we choose the normalization a = 1): 



H 2 (z) = 



2Vp 
6- A 2 



b 2 {z + l)- (aX+3 ^ +a 



(24) 



In redshift variable the energy density of DM is p m (z) = 
(M/X )(z + l)« A + 3 7 = p ( z + i)"A+3 7j so the DM di _ 

mensionless density parameter can be given as function 
of z also: 



n m (z) 



(6 - A 2 )p (z + l)"A+3 7 -A 2 
2V b 2 (z + l)"A+37-A 2 + 1/52 



(25) 



In order to constrain the parameter space of the so- 
lution one should consider the model to fit the obser- 
vational evidence about a universe with an early mat- 
ter dominated period and a former transition to dark 
energy dominance. Therefore A 2 — aX — 3j < 0, be- 
sides, one can write Q m (oo) = (1 — e), where e is a 
small number (the small fraction of dark energy compo- 
nent during nucleosynthesis epoch). In correspondence 

b 2 = 2 6 (7-e)Vh • Taking into account the observational 
fact that, at present, (z = 0); £! m (0) = 1/3, then 
1/b 2 = 2 — 3e. After this Eqn.(25) can be rewritten in 
the following way: 



Q m (z) = (l-e) 



+ 



aA+3 7 -A 2 



(z + l)«A+3 7 -A2 + ( 2 _ 3e ) 



(26) 



Other physical magnitudes of observational interest are 
the quintessence equation of state (EOS) parameter and 
the deceleration parameter 



W0 



-1 + 



A 2 /3 
l-O, 



respectively. 



-1 + — + ^O ro , 



2. Case B 



(27) 



(28) 



A second very simple choice is to consider an exponen- 
tial coupling function of the form: 



A(2 — ,) 

X = X e^^< 



In this case equation (12) simplifies to 
dV 



XV = 0, 



(29) 



(30) 



which can be easily integrated to yield to a single expo- 
nential potential: 
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V = Voe-^. (31) 
In consequence the equation (15) can be written as 



/ 



dw 



Vw 2 + A 2 



fJ.(T + To) 



(32) 



where w = cxpf- 7 ^/^ ), A 2 = V X /M and fi = 
r jy/M{\ 2 — 6)/8, so we have the explicit solution 

0(t) =0 o + ln(sinh[^(r + r o )]^ra) , (33) 

and, consequently: 

4A 

o(t) = a sinh [/x(r + r )] tC 6 -^) . (34) 

The dimensionless density parameter and the Hubble 
expansion parameter are given, as functions of the red- 
shift, through: 



fi m (z) 



6 - A 2 (1 + z )37+7(6-A 2 )/2-A 2 
~6A^~ (1 + z )3 7 +7(6-A 2 )/2-A2 + A 2 ■ 



(35) 



and 



H(Z) = B\j^(l + z )37+7(6-A 2 )/2 + (1 + z) X^ ( 36) 



respectively, where _B = ^/2Uo/(6 — A 2 ). Note that, O m 
is a maximum = (6 — A 2 )/6A 2 ) at z = oo. Taking 
into account that observations are not exact we admit 
that, at the epoch of nucleosynthesis £l m (oo) = (1 — e), 
where, as before, e is a very small number (e = [6(A 2 — 
1) + A 2 ]/6A 2 ). Then, if we assume (as before) that, at 
present, fi m (0) = 1/3, then (35) can be rewritten as 



£l m (z) = 



(1- £ )(1 + z) 3 ^ 6 - a2 )/ 2 - a2 

(1 + z )37+7(6-A2)/2-A^ + ( 2 _ 3e ) 



(37) 



Other physical magnitudes of observational interest arc 
the quintessence EOS parameter and the deceleration pa- 
rameter, that are given in equations (27) and (28), re- 
spectively 



IV. EXISTENCE AND STABILITY OF THE 
SOLUTIONS 

In this section we do not aim at a complete study of the 
stability in the case of the double exponential potential. 
Instead, we use the following simplifying fact. Note that, 
in terms of the scale factor, the potential (17) can be 
written as follows: 



V(a) = V a- x2 + W a"( aA+37 \ 



(38) 



consequently, at early time (small a) the second term in 
the right hand side (RHS) of (38) prevails so, the first 



term can be neglected and we can talk about a single ex- 
ponential potential with the constant — (a + 37/A) in the 
exponent (see Eqn.(17)). At late time (a 3> 1), instead, 
the first term in the RHS of (38) dominates, so we can 
talk about a single exponential with the constant —A in 
the exponent. For this reason, and owing to the fact that 
in both cases (case A and case B), the coupling function 
is of an exponential form, we focuss (as first approxima- 
tion) on the stability analysis of the model with a single 
exponential potential of the form: 

V = Vbexp(-at^), (39) 

and the following coupling function: 

X = X Q cxp(60). (40) 

The constants a and b depend on the model. Using the 
variables [9] 



V 



V6H 



, y : - 



V3H 



(41) 



the system of dynamical equations (3-5) can be arranged 
in the form of an autonomous dynamical system: 



' - 2 2 2\ 

x =a \J 2 V ~ b \J 2*- V > 

-3s + \x\2x 2 + 7(1 - x 2 - y 2 )] 



V 



| 2 /[2a; 2 +7(l 



y 2 )} 



(42) 



<'\l^yx (43) 



where, as before, the prime denotes derivative in respect 
to the variable N = ln(a). If in equations (42) and (43) 
we set b — 0, the case studied in Ref. [9] is recovered. 

Using the new variables, the Friedmann equation looks 
like: 



Pn 



3H 2 



(1 



y 2 l 



otherwise: 



2 1 2 

x +y . 



(44) 



(45) 



Since p m > 0, the new variables are confined to the 
compact region (see [9]), < x 2 +y 2 < 1 so, the evolution 
of the system under investigation is completely described 
by trajectories within the unit disc. The lower half-disc, 
y < 0, corresponds to contracting universes, but we are 
interest in expanding universes. On the other hand, due 
to the following discrete symmetries : (x,y) — ► (x,—y) 
and t — > —t, we consider the upper half-disc, y > in the 
subsequent discussion [9]. 

The effective EOS parameter for the quintessence field 
at any point is given by 



70 



2x 2 



(46) 



P4> V{(j))+<t> 2 l2 x 2 +y z 

We are able to find up to five fixed points (critical 
points) where x' = and y' = 0, which are listed in 
Table I. 
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X 


y 


Existence 


Stability 


&<4> 




-1 





All a and fo 


Stable node for a < — \/6 and 

b < y/l ( 7 - 2) 

Unstable node for a < — \/E and 

b< v/f( 7 -2) ^ 
Saddle point otherwise 


1 


2 


1 





All a and fo 


Stable node for a > \/6 and 
5 > 71(2-7) 
Unstable node lor a < yo and 

6< 7f(2-7) 
Saddle point otherwise 


1 


2 




n 
U 


A 11 ^ T, ^ . / 3(7~ 2) 


c ^ u 1 ^ ^ ^ t nv . 7; ^ 37(2—7) 


2b 2 


O 

z 


7-2 


All a ana o <i \ — 2 — 


otaoie node tor a < — ^ — 
Saddle point for a > 


3(7-2) 2 


a 

s/6 


\A-f 

V 6 


All a and fe 


Stable node for a 2 < 6 and 6 > " 2 r 37 

Unstable node for a 2 > 6 and 6 < " T 37 
Saddle point otherwise 


1 


3 


fS 7 


/l "1 3 7 3 7 2 


All a and 6 < "'l 37 


Stable node or Saddle point 
(see analysis in table (II)) 


37+S 2 -aS 


37 2 


V 2 a-i) 


Y a-S> ~ (S-6) 2 2(a-b) 1! 


— a 


(a-b) 2 


37 + b 2 -ab 



TABLE I. Existence and stability of the critical points. 



Case 1 


0<a<\/6A0<5<aA0<7<l-7i 


Case 2 


a = \/6A{0<&<&iA0<7<l-7iVbi<&<aA0<7<72} 


Case 3 


\/6 < a < 2 v / 2 A {0 < 6 < & A (0 < 7 < 1 - 71 V 1 + 71 < 7 < 2) V 60 < 6 < 61 A (0 < 7 < 1 - 71 V 1 + 71 < 7 < 72) 
V6 = 61 A (0 < 7 < 1 - 71 V 1 - 71 < 7 < 72) V 61 < b < a A < 7 < 72} 


Case 4 


a = 2 v / 2 A {0 < 6 < 60 A (0 < 7 < 1 - 71 V 1 + 71 < 7 < 2)V 
6=&oA(0<7<l-7iVl-7i<7<2)V&o<&<aAO<7<72} 


Case 5 


a > 2V2 A{0<6<6iA(0<7<1-7iV1+7i<7<2)v6 = 6iA(0<7<1- 71V 
l-7i<7<2)Vbi<6<6 u A0<7<2v6()<b<aA0<7<72} 



TABLE II. Cases in which the point 5 is a saddle point (otherwise it is an stable point). We have used 71 = — 
72 = § (a 2 — ab), bi = § — + a 2 and 60 = ~ 6 ^ a ■ A and V are the logical operators "and", "or" respectively. 



A. Double exponential potential 

In this case the constant b = a is fixed. Since we have 
used the constrain (9), then the dynamical variable x is 
fixed to be x = ^= so, depending on a, we are led with 

two possibilities: i) a = A; the point 4 in table I is an 
attractor quintessence dominated regime (Fig.l) and ii) 
a = a + 37/ A; the point 5 in table I is a matter-scaling 
attractor (Fig. 2). Therefore, we can trace the evolution 
of the universe from a matter-scaling attractor in the past 
into a late time, DE dominated attractor. 



B. Single exponential potential 

In this case two critical points emerge. Actually, al- 
though a and b are fixed (we have a single exponential in 
this case): a = A and b = A/2 (2 — 7), the points 3 and 
4 in table I are consistent with the constrain (9). The 
point 3 is a scaling solution and the point 4 is a solution 
dominated by the scalar field. 

As in the former case, the evolution of the universe 
transits from a scaling attractor into a scalar field domi- 
nated attractor in the future. This fact makes possible to 
deal with the fine tunning and the coincidence problems 
[5,6]. 
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FIG. 1. The phase plane for 7 = 1, A = 0.245 and a = 1. 
The point (0.1, 0.99) is an attractor. All trajectories in phase 
plane diverge from the unstable nodes at (1, 0) and (-1, 0). 
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x 

FIG. 2. The phase plane for 7 = 1, A = 0.245 and a = 1. 
The point (0.81, 0) is an attractor. All trajectories in phase 
plane diverge from an unstable node (-1, 0). 

V. OBSERVATIONAL TESTING 

Although we can not talk literally about observational 
testing of the solutions found (we point out that the ob- 
servational facts are not yet conclusive in some cases and 
are not accurate enough in others), in this section we 
shall study whether our solutions agree with some obser- 
vational constrains that are more or less well established. 
The main observational facts we consider arc the follow- 
ing [13]: 

1. - At present (z = 0) the expansion is accelerated 
(9(0) < 0). 

2. - The accelerated expansion is a relatively recent phe- 
nomenon. Observations point to a decelerated phase of 
the cosmic evolution at redshift z = 1.7. There is agree- 
ment in that transition from decelerated into accelerated 
expansion occurred atazw 0.5 [14]. 

3. - The equation of state for the scalar field at present 
uj,p(w v ) ~ —1 (it behaves like a cosmological constant). 



With a 95% confidence limit v < -0.6 [15]. t 

4. - Although, at present, both the scalar (quintessence) 
field and the ordinary matter have similar contributions 
in the energy content of the universe (Q m (0) = 1/3 =^> 
^^(O) = 2/3), in the past, the ordinary matter domi- 
nated the cosmic evolution, * meanwhile, in the future, 
the quintessence field will dominate (it already domi- 
nates) and will, consequently, determine the destiny of 
the cosmic evolution. 

5. - As stated before [18], nucleosynthesis predictions 
claims that at 95% confidence level. Q v {lMeV ~ z — 
10 10 ) < 0.045. 

6. - During galaxy formation epoch [19], around z ~ 
2 — 4, the value of quintessence density parameter is Q v < 
0.5. 

Now we proceed to " observationally" test the solutions 
found in the cases studied in the former sections. 



A. Double exponential potential 

In this case the solution depends on two parameters, 
A and e (a and 7 are fixed: 7=1, meaning cold dark 
matter dominance at present, and a = 1). To constrain 
the values of these parameters we use the aforementioned 
facts (1-6). Using these constrains we bound the parame- 
ter space (A and e) by means of a computing code. These 
parameters range as < e < 0.045 and < A < 0.38. 
Then we select a pair from the above interval in order 
to compare our model with the supernova data. We fol- 
low the procedure used in [20] to further constrain the 
parameters in the model. § For instance, if we consider 
e = 0.01, the x 2 distribution has a minimum in the to 
direction** at too = 24; however, it has no minimum in 
the A-direction, meaning that the parameter A can take 
any values in the interval < A < 0.38. In the Figure 
3, x 2 is plotted as a function of the free parameters A 
and Too (we chose e = 0.01, A could be any value in the 
physically meaningful range). Although the model fulfils 
the observational requirements (1-6) given at the begin- 
ning of this section, we point out that other observations 
should be considered in order to further constraint the 
space of parameter. 



^Corasaniti and Copeland [17] found that the determination 
of the third peak in the BOOMERANG data limits the value 
today of equation of state — 1 < ui v (0) < —0.93. 

*A sufficiently long matter dominated decelerated phase is 
needed for the observed structure to develop from the density 
inhomogeneities [16] 

§ In fact this procedure cannot be used to further constrain 
the space of parameter, since the free parameters in our model 
are not very sensitive to these observations. 

"mo is a parameter connected to the absolute magnitude 
and the Hubble parameter. 
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0.5 22 

FIG. 3. The function \ 2 (C in the fi gure) is plotted vs 
the free parameters m and A (L in the figure), for the model 
with double exponential potential. We chose the parameter 
e — 0.01. Note that A could be any value in the physically 
meaningful range so, SNIa luminosity observations do not al- 
low for further constrain of the parameter space. Other ob- 
servations could be considered for this purpose. 

B. Single exponential potential 

In this case our solution depends on two parameters 
also A and e (recall that we fix 7 to be unity). We apply 
exactly the same procedure as that in the former subsec- 
tion to constraint the parameter space. We found that 
the physically meaningful region in the parameter space 
is bounded by < e < 0.045 and < A 2 < 0.135. 

We want to stress that, as in the former case, the free 
parameters of the model are not very sensitive to the 
SNIa data, so, the standard used, for instance in [20], is 
not suitable for further constraining the space of param- 
eters of the model. Other observational evidence could 
be considered for this purpose. 

It should be pointed out that, in two cases (A and B), 
the meaningful region in parameter space is chosen such 
that the main observational facts (1-6) explained at the 
beginning of the section are fulfilled. As an illustration, 
in Figure 4, we show the evolution of both dimension- 
less DM and scalar-field energy densities Q m and Q^, re- 
spectively vs z for the single exponential potential (the 
Case A is very similar). In Figure 5, the evolution of the 
equation of state is shown vs z, while in Figure 6, we 
plot the deceleration parameter to show the transition 
redshift when the expansion turns from decelerated into 
accelerated (z s» 0.55). 

VI. CONCLUSIONS 

We have found a new parametric class of exact cos- 
mological scaling solutions in a theory with general non- 
minimal coupling between the components of the cosmic 
mixture: the cold dark matter and the dark energy (the 
quintessence field). To specify the general form of the 
coupling we were inspired in a scalar tensor theory of 



redshift 

FIG. 4. The evolution of Q m (thick solid line) and 
(dashed line) vs z is shown for the model with a single expo- 
nential potential. The following values of the free parameters 
e = 0.01 and A = 0.3 have been chosen. Equality of mat- 
ter and quintessence energy density occurs approximately at 
0.3 -0.4. 

1 l I 
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FIG. 5. We plot the dynamical EOS parameter of the scalar 
field vs z for the Model B (single exponential potential) . The 
values of the free parameters chosen are e = 0.01 and A = 0.3 
respectively. Note that the DE evolves from almost being 
attractive matter (dust in particular) in the past to behaving 
as an almost a cosmological constant at present. 

gravity written in the Einstein frame. We have studied 
particular coupling functions that lead to self-interaction 
potentials of the following class: i) double exponential 
potential, ii) single exponential potential. In order to 
derive exact flat FRW solutions we have assumed a lin- 
ear relationship between the Hubble expansion parame- 
ter and the time derivative of the scalar field. The sta- 
bility and existence of these solutions have been studied. 

In Case A (double exponential potential) the universe 
evolves from an attractor regime, characterized by the 
scaling of matter and dark energy, into a late time, dark 
energy dominated attractor regime. In Case B we see 
that a single exponential can lead to one of the follow- 
ing scaling attractor solutions [21]: either 1) a late-time 
attractor is a state when the scalar field mimics the evo- 
lution of the background (DM) fluid, or 2) the late-time 
attractor is the scalar field dominated solution. We show 
that the assumed linear relationship between the Hub- 
ble parameter and the time derivative of the scalar field 
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FIG. 6. We plot the evolution of the deceleration parameter 
q vs 2 for the Model B (single exponential potential). As be- 
fore, the following values of the free parameters e = 0.01 and 
A = 0.3 have been chosen. The transition redshift when the 
expansion turns from decelerated into accelerated occurred at 
z T w 0.54- 0.56. 

always leads to attractor scaling solutions. 

In all cases, the models were tested against the widely 
accepted observational facts. The space of parameters is 
then reduced to a compact region whenever it is possible. 
It is shown that the existing observational evidence com- 
ing from SNIa data is not enough to further constrain 
the space of parameters. In all cases the models fit with 
enough accuracy the observational data. 

We conclude that models with non-minimal coupling 
between the dark energy and the dark matter are easy to 
handle mathematically if one assumes a suitable dynam- 
ics. Since, in the cases studied, the coupling function is 
an exponential (these differ only in the constant in the 
exponent), we point out that the resulting theory belongs 
in the class of Brans-Dicke-type of theory, written in the 
Einstein frame. 
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